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In this paper we study the relation between the light-front (infinite momentum) and rest-frame 
descriptions of quarkonia. While the former is more convenient for high-energy production, the latter 
is usually used for the evaluation of charmonium properties. In particular, we discuss the dynamics 
of a relativistically moving system with nonrelativistic internal motion and give relations between 
rest frame and light-front potentials used for the description of quarkonium states. We consider 
two approximations, first the small coupling regime, and next the nonperturbative small binding 
energy approximation. In both cases we get consistent results. Our results could be relevant for the 
description of final state interactions in a wide class of processes, including quarkonium production 
on nuclei and plasma. Moreover, they can be extended to the description of final state interactions 
in the production of weakly bound systems, such as for example the deuteron. 


I. INTRODUCTION 


In the general case, a direct boosting of the wave function components from one frame to another presents a 
complicated dynamical problem, which mixes states with various parton content |l], However, for certain systems, 
for example positronium, heavy quarks systems, or a nucleus with vanishingly small binding energy, the Fock state is 
dominated by the lowest component with minimal number of partons. The dynamics of such systems is described in 
potential models, and the wave functions in different reference frames can be related to each other: it is well-known 
that in the light front quantization approach the wave function is invariant with respect to longitudinal boosts M, 
so a transformation of the rest frame wave function requires just to rewrite it in terms of light-front variables 4|. 
This solves a boosting problem for a class of processes, in which due to factorization theorems only the light front 
wave function of the initial or final hadron is needed. However, in certain cases a more sophisticated approach is 
necessary: when, for example, processes inside extended objects like nuclei 11 and plasmas [7jj 1 are considered. In 
the evolution of such systems, an effective potential should be taken into account at all stages. In the light front, which 
is a natural choice for description of high energy processes, usually the quark potentials are much less understood 
than in the rest frame. For example, nothing is known about the light-front potentials at nonzero temperatures. The 
main goal of this manuscript is to fill this gap and provide a method which would bridge the light front and rest frame 
approaches. For this purpose we use the relation between the light front and the so called infinite momentum frame. 

In what follows, for the sake of definiteness we will consider the particular case of quarkonium (charmonium 
or bottomonium), tacitly assuming that results could be extended to other systems. Quarkonium has been very 
well understood in the rest frame |21| . even for nonzero temperature [22b 29). Experimentally there are data for 
quarkonia production both on the proton and nuclear targets (see e.g. f2oT~ 3oT~311 for review and references therein). 


In photo- and electroproduction on protons, factorization theorems hold, and as was mentioned above, the final state 
wave function can be obtained jj. In production on nuclear targets, as was discussed in [fr, d], the contribution of 
final state interactions (FSI) is important and requires knowledge of the light-front potentials. Recently the relation 
between the rest frame and light front quarkonium potentials has been studied in AdS/QCD framework f32l . l33j . 
Basing on equality of the spectra, it was suggested that the transformation could be nonlinear, albeit inclusion of the 
higher order 0(l/m q ) terms inevitably has some ambiguity and differs between different authors 


From our 

point of view, inclusion of such terms is not justified since at the same order there are contributions from the omitted 
multiparton Fock states. 

In this paper we re-visit the problem of bound state of two heavy quarks, where the whole system is moving with 
relativistic momentum in the laboratory frame. In Section [H] we analyze the dynamics of the system perturbatively in 
the moving frame and reduce a Bethe-Salpeter equation to a Schroedinger equation (0) for a nonrelativistic internal 
motion of a moving dipole. The latter is preferable for many practical applications because the near-onshellness 
of both heavy quarks in a Bethe-Salpeter equation complicates the numerical treatment [34]]. The novelty of our 
approach is that, in contrast to existing treatments, from the very beginning we consider a moving quarkonium. Also, 
we give a generalization of the pNRQCD lagrangian © in a moving frame. 

In view of the fact that for charmonium a s ( m c ) « 0.25 is not very small, in Section IlIIl we develop a different 
approach to the non-relativistic system, based on analyticity and unitarity constraints, in which we assume smallness 
of the parameter e/A qcd and which is not based (at least explicitly) on the smallness of the QCD coupling. In this 
approach, we get exactly the same Schroedinger equation (JHl). In this method of derivation we pursue a practical 
aim: to formulate a set of rules on how to use the phenomenology developed for the description of J/'F-meson in the 
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Figure 1: Bethe-Salpeter kernel for the two-parton system. 


laboratory frame, in any moving reference frame. 


II. SCHROEDINGER EQUATION FOR HEAVY QUARKONIUM IN A MOVING REFERENCE 
FRAME, FROM THE BETHE-SALPETER EQUATION 


As it is well known @,53, in the rest frame there is a natural small parameter, the velocity of the heavy quark 
v -C 1. A systematic expansion over this parameter leads to an effective theory, NRQCD. In a perturbative pNRQCD, 
(which requires an additional assumption a s (m q ) <C 1) the velocity of internal motion of the bound quark v ~ a s , so 
a small-'U expansion is equivalent to a systematic perturbative expansion. For asymptotically heavy quarks, there is 
a hierarchy of scales 

m q > m q a s > m q a 2 s » Aq C d (1) 


where m q is the mass of the heavy (charm or bottom) quark. The scales in this hierarchy have a straightforward 
physical meaning: the soft scale m g o. corresponds to inverse size of the quarkonium system; the ultrasoft scale m q a 2 
is the typical binding energy e, etc @,53- 

The dynamics of the system is described by the Bethe-Salpeter equation (BSE), which is explicitly invariant in any 
system and contains both a rest frame and a light front Schroedinger equations as special limits. In an explicit form, 
the BSE is written as (see Fig. [T]) 


iXab (p~,-p+) = S ac ( P -) 


d 4 q 

(2A 


Acejd{p-,q-p+-p- 1 q-p+) Xef ( q , q-p- -p+) S d b(-p +), 


( 2 ) 


where p _ and — p + are momenta of the quark and antiquark respectively. In Eq. Q A ce jcL denotes the 2-quark 
scattering amplitude. Recall that in this equation all ingredients are relativistic invariants and do not depend on the 
reference frame. 

In a perturbative QCD, the propagator S in the leading order should be replaced by the free quark propagator, 
and the vertex part A should be replaced by a single-gluon exchange in the f-channel, so fS|) simplifies to 


ix{P-,-P+) = S{p-) 


A (q, q-p-~p+)iv o ( , 

A) 4 (g — P-) 2 + *0 1 P+> ’ 


( 3 ) 


where from now on we omit explicit Dirac indices. If we denote po = P/2, where P is the momentum of the 
quarkonium, and define small deviations as Sp- = p- — po, —5p+ = — p + — po, after some well-known algebraic 
manipulations (see details in Appendix [Aj , we may reduce (J3J to 


ip-,-p+) 


g{5p-)g{-8p+) 


/ 


A 


x(g, q-p- -p+) 

(q-p -) 2 + i o 


X = a +7mX7A- 


( 4 ) 

( 5 ) 


Note that while \ is a matrix, equation Q is a scalar equation (doesn’t mix components), and for this reason in what 
follows we may omit the spin structure of \ and treat it as a scalar function. Taking into account that the projection 
of the vector q—p on P is O (a 2 ') (this corresponds to a generalization of the rest frame instantaneous approximation 
k 0 ~ 0), after some algebra (see Appendix 0 we can obtain a Schroedinger equation for the internal motion in the 
form 




^ (C, r_ l) = V (C,r ± ) V (C, rj_) 


( 6 ) 
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where we defined a wave function df as \II = f dX x (q, 9 — P)^_ p , a potential 1/ is proportional to the trace of the 
gluon propagator, V = g^IL^/A, where n M „ is the gluon propagator, and the derivatives are defined as 

d ,f = — (cosh r/do + sinh 77 ^ 3 ), (7) 

= — (sinh 77^0 + cosh 77 ^ 3 ). ( 8 ) 

Here we introduced the shorthand notations 77 = In (P + /M), and M for the mass of quarkonium. The relation of the 
effective potential V to the rest frame potential will be discussed in the next section. Equation © provides a smooth 
interpolation between the rest frame (77 = 0) and light-front frame (77 —>• 00 ). This equation can be also obtained from 
the pNRQCD lagrangian 0 , which in a moving frame gets the form 


L 


= -4 (ids + — ~ — - Vs) + 4 (iD( + — ~ — - V 0 \ + 

y m q m q J u y m q m q J 

+ -^-V A (r) (4<P ^^ 0 + + 

+ 9 ^- (4 <p^^ 0 + 4 <p"p^o) -Jp^p;,, 


( 9 ) 


where H /5 and H/q = A>Qt a are the fields of the singlet and octet qq-p&ir (T in © corresponds to in ©), is 
the momentum of the moving quarkonium, + ig[A^, \I>o], A M is the gauge field, and Vs(r), Vo{s), Va{t) and 

Vb(t) are the potentials for singlet, octet diquarks and a transition matrix elements. While in pQCD they are given 
by well-known perturbative expressions, for a quarkonium propagating inside matter they become more complicated 
and get an absorptive part. A detailed study of the absorption mechanism is out of scope of this paper and will be 
presented elsewhere. 


III. SCHROEDINGER EQUATION FROM ANALYTICITY AND UNITARITY 


Notice that the single-gluon exchange in f-channel gives only a Coulomb term in the effective potential. As we 
have discussed in the introduction, such an approach can be justified only for very heavy quark states (say bottomo- 
nium), whereas for the JJJS meson higher order corrections are essential. This fact manifests itself in the developed 
phenomenology (Tll - ITsl . [23l . [24j in the rest frame, in which an additional confining potential is added to one gluon 
exchange. This extra term is generated by nonperturbative interactions and introduced into the Bethe-Salpeter either 
as additional vector or scalar f-channel contribution (see e.g. El). 

Bearing this experience in mind, in this section we are going to obtain the Schroedinger equation for the non- 
relativistic system without using the explicit form of A in Eq. J2J . In the following discussion, it is convenient to write 
the amplitude % as a function of the variables s = (p_ — p + ) and p = p_ + p+ instead of parton momenta p± . The 
amplitude x is analytic outside the real axis, so it can be represented as 



Imx(s', p) 
s' -M 2 + i0' 


( 10 ) 


where M is the mass of the qq meson. The imaginary part of the amplitude x can be evaluated directly from 
using the unitarity constraints 


d 3 q x(s', q) A 2 ^ 2 (s', ( q - pfj 6 (s' - (2 q - p) 2 ) + ^ Xi-m &) A n ^ 2 

n =3 


(ii) 


In this equation A 2 _ ¥2 (s', p — p') is the 2-particle scattering amplitude. If we neglect the contributions of higher Fock 
states (n > 3 in Eq. (fTTT) 1 in the large- m q limit, we can re-write Eq. (fTTfll and Eq. (HID in the following form 


* ^ Ma ) = h / ^ L *' S-M’ + iO (»'.<«- P) 2 ) * (»' - < 2 « - tf) + 0 4 ) . ( 12 > 


where ma ~ 1 GeV is the effective mass of gluon E3- These estimates stem from the assumption that the second 
term in Eq. m has the same order of the magnitude as the first one, but its contribution starts to be essential for 
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s' > (2 m q + me) 2 - In the dispersion integral the first term shows the enhancement oc 1 /(Ale), while the second has 
the denominator ex 1 /(M m g ). It is difficult to make better estimates for accuracy of our calculation due to the lack 
of understanding of non-perturbative QCD. 

It should be stressed that for very heavy quark-antiquark systems, at small distances the omitted terms have an 
additional suppression of the order of 0(a 2 (m q )). In the case of realistic system such as the J/'k-meson we cannot use 
this smallness, and in order to evaluate the highly excited states in Eq. m we have to rely on models. It turns out 
that in phenomenological models for the J/’I'-meson, a substantial contribution stems from the string-like potential. 
Indeed, V e g(rj/^) — IpQCD (Aj/vO ~ cr (o/V’) ~ 0.4 GeV for the string tension er = 1 GeV/fm. Bearing this fact in 
mind we can estimate the contribution of the multi-particle states in the unitarity constraint, by comparing the mass 
of the next string excitation or, in other words, the mass of the next resonance on the J/’k Reggeon trajectory with 
mass Mr. This contribution is suppressed by the parameter 


Q 


2 Me 

M 2 r - M 2 ' 


(13) 


The binding energy e is a poorly defined object, since the mass of the heavy quark is scheme-dependent. For 
the case of charm quark, the estimates for the mass m c vary between 1.27 GeV and 1.8 GeV, so as an upper value 
estimate, we take e max ~ 2Mr> — Mj/^p ta 600 MeV. Taking M = Mj/^p and Mr = 3 686 ) [3y| for the first excited 

state in a channel with J/iJ) quantum numbers, we get an upper limit for the parameter Q max ~ 1/2. In the small-e 
approximation we may simplify (1121) as 


A 2 ^ 2 (s', (q - p) 2 ) « A 2 -> 2 (m 2 , (q - p) 2 ) • 


(14) 


It should be mentioned that the assumption of Eq. (070) means that amplitude A 2 -> 2 (s', (q — p)j has no singularities 
related to the quark-antiquark state, and this amplitude can be viewed as a sum of quark-antiquark irreducible 
Feynman diagrams in QCD. Introducing a new function 


tf(s, q) = 


X'(s, q) 

s — M 2 + iO ’ 


(15) 


we may rewrite (1121) in the form 


(p 2 - M 2 ) 4/(M 2 , P) = ^J d 3 qV(q-p)^(M 2 , p). 


(16) 


In the instant form we may split the vector of relative motion p into a part p| | collinear to the quarkonium momentum 
P and part Sp orthogonal to it. After some trivial algebra, (fTHl) reduces to a simple Schroedinger equation ©. In the 
front form (1161) reduces to the result obtained in fH, 


MH 


( Pj_ + ITlq 
y x(l — x) 





(17) 


where M is the mass of the meson, q± is the quark transverse momentum, x is the light-front fraction of the longitudinal 
momentum of the quark, and U is the operator of the potential energy. In the heavy quark mass limit, there are 
two natural small parameters, O (a s (m q )) and O (Aqcd/w 9 ). Typical momenta of the quarks in a nonrelativistic 
system are of order m q O (a s ( m q )), whereas the binding energy is e ~ m q O (a 2 ( m q )). Taking into account that the 
light-front fraction x ~ 1/2 + Sx 1 where 5x s=s O (a s ( m q )), we can see that (TT71) reduces to 

eU/ = f -p 4 to § x 2 _|_ -vp (18) 

\m q Am q J 

If we make a Fourier transformation of 4' to coordinate space according to 

rx) = J (* = \ +<y®, gx) (19) 


1 


Recall that fraction x is invariant with respect to boost in the longitudinal direction and can be calculated in the rest frame, x = 
(<70 + q z )/M = ^ + (e + y / mge)/M, where M is the mass of the quarkonium. 
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where £ is a parton separation along the z -axis, and <f> satisfies the ordinary Schroedinger equation 


e</>(C, Cl) = 


Au - d 2 , . \ 

-4VU((, rj. , 

m q J 


( 20 ) 


where V = U/ (4 m q ). We would like to emphasize that in contrast to |33j the potential V is local only in the coordinate 
space; in a LF frame the interaction has a form of a convolution 


V</>(x, r±) = V (x-xi, r±)(j){xi, r±)dx\. 


( 21 ) 


The Green function of the internal motion in the infinite momentum frame satisfies an evolution equation 

G(x + ,C, r±;xf,Curi,±) = 


M dx + 
-A x - d? + 


( 22 ) 


L G(x + , £, r±;xf,(i,ri,±) +VG (x + ,£, r±\xf ,Ci,n,x) , 


Finally, we would like to address how the light front and rest frame potentials are related to each other. From 
Eq. m we see that this amplitude is a function of only four-dimensional transfered momentum. As we can see 
from mm , a relation between the potentials in different frames can be extracted from the transformation of the 
gluon momentum k = q—p . In a rest frame, the vector k has a negligible ko component (instantaneous approximation), 
so k 2 Rj —k 2 + O (af) . In an infinite momentum frame, we may rewrite 


k 2 = ( q-pf = (x- xi) 


m„ 


■pi 


m„ 


■9x 


Xi 


d. 




(23) 


where light-front fractions are defined as x = p + /P + and X\ = q + /P + . Note that x — aq ~ O (qj_/m q ), so both 
terms are of the same order of smallness and the omission of the terms in a previous line is justified due to extra 
~ 0(q±/m q ) suppression. After some lengthy but straightforward evaluations (see details in Appendix [Bj we may 
get for the interaction term in the infinite momentum frame 

(x, r±) = J dxiK (x, xi, r±) 4' (x lt r ±), (24) 

where the kernel 


K (x,xi,rj_) 



exp (2 im q (x — x±) z). 


(25) 


In Figure [2] we compare the ground state (J/VO wave function evaluated from (1221 ) with the one obtained with the 
phenomenological prescription |4j. For the sake of definiteness, in both cases we used a rest frame potential from [28l . 
As we can see from the plots, for the symmetric configuration (x « 0.5) the wave function given by a prescription [J] 
coincides within the accuracy of numerical evaluations with the eigenfunction of (1221) . However, for large-a; there 
is some disagreement between the two wave functions. This could be easily understood, taking into account that 
x — 1/2 is a small O (a s ) parameter which is neglected. The average width of the charmonium wave function is 
(l x — ll) ~ (ps/2 rn q ) as 0.2, which signals that the relativistic effects, albeit small, are not completely negligible. 


IV. CONCLUSIONS 

In this paper we re-visit the problem of the wave function for the heavy quark-antiquark system in an arbitrary 
reference frame. In two regions: (i) m q as (m q ) A.qcd and (ii) e/Aqcd 1, we constructed a Schroedinger 
equation for the internal motion of J/ip moving in the laboratory frame with relativistic center-of-mass momentum 
(see Eq. (O). In these cases, this equation reduces to the rest frame and infinite momentum Schroedinger equations, 
and thus gives a smooth interpolation between the two limits. We studied a relation between the rest frame and 
light-front potentials of quarkonium. We found that the two are related by a linear transformation (1241) . which 
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Figure 2: Comparison of the ground state wave functions evaluated with (1221) (solid red line) and with a prescription from [3] 
(dashed blue line). Upper row: r^-dependence at fixed x. Lower row: x-dependence at fixed values of the parameter rj_. 


corresponds to a large-m 9 limit of a nonlinear transformation in [32l . [3dl | . The omitted terms are of higher order 
in 0(Aqcd j (cts (nt g )) and O (e j A qcd'j and their inclusion is not justified in the two-parton approximation. We 

checked the phenomenological prescription Q and found that in the region x ss 0.5, which is relevant for charmonium- 
related problems, it is well justified. Also, we got that in the infinite momentum frame a potential is given by a 
convolution (1241251) in light-front variable x, in contrast to what was found in [33]. 

Our results could be relevant for the description of the final state interactions of a wide class of processes includ¬ 
ing quarkonium production on nuclei or plasma. Besides, they could be extended to the description of final state 
interactions in the production of weakly bound systems, such as for example the deuteron. 


Acknowledgments 


This work was supported in part by Fondecyt (Chile) grants No. 1130543, 1140390, 1140842 and 1140377. We are 
grateful to Valery E. Lyubovitskij for discussion of potential transformation properties in the AdS/QCD framework 
and for providing the reference [33]. 


Appendix A: Derivation of the Schroedinger equation 


The Bethe-Salpeter equation (BSE) has a form 

d 4 q 


iXab{p-,~P+) = S a 


(p-)/ 


(2tt)< 


A C ejd (p~,q — p+ -p-,q-p+)xef {q, q-p- - p+) S d b (-??+), 


(Al) 


where p_ and — p + are the momenta of the quark and antiquark respectively, and we explicitly have shown lower 
Dirac/flavour indices. The Bethe-Salpeter amplitude defined as 


Xab 


(P 1 ,P 2 )= J d A x 1 d 4 x 2 e i ^- p ^(0\^ a ( Xl )j b (x 2 )\P) 
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is not a wave function, but can be projected to a wave function integrating along some time-like curve. The propagator 
S in leading order may be replaced by the free quark propagator, and the amplitude A is dominated by a single-gluon 
exchange in the f-channel. Note however that single-gluon exchange produces just a Coulomb potential, whereas in 
realistic models we have also a string, so we assume that the gluon propagator has a form = g^V, where V is a 
Fourier image of the potential. In this case we may rewrite the equation (IA1I) as 


ix(P-,~P+) = s 



d 4 q 


V{q- p -) 7 VX (<?, Q-P-- P+) 1»S (-P+), 


(A2) 


Starting from this equation we omit the Dirac indices since we have a matrix identity. If we denote po = Pj/^/2 
and define small deviations as Sp- = p- — po, —Sp+ = —p+ — po, we may expect that a deviation of quark momenta 
from po are small; the terms ~ Spi may be neglected in the numerator but should be retained in the denominator, so 
we get 


S(p-) = 


P- + m q po +m q 


P- - ra, + *0 2 m q po-Sp- _ , -q 

m a 2m a ' 


= A +g (Sp-) 


(A3) 


C ( ^ - p+ + rn q po-m q 1 . , , s 

s (~P+) = -5—^ ~ -75- - -7- 73 -= A ~g - Sp+ 

2m n P 0 -Sp+ _ °p+ 


P+ — rriq + zO *m q po-qp+ "r+ | 


(A4) 


m n 2m„ 


where we introduced a shorthand notation for the projectors 

Po ± rn q 


A± = 


2 m a 


and a new function 


9(Sp) 


Po-Sp 


Sp 2 
2m q 


•i0 


Multiplying both parts of (I A2h by A+y^ from the left and by y„A_ from the right, after some manipulations with 
Dirac algebra, we get 

/ j4 

—^v(q-p-)x(q,q-p--p+), (As) 

X = A+y^yy^A. (A6) 

In a moving system, a projection of the gluon momentum q — p onto the direction of vector po is O (o^), whereas 

all the other components are O (a s ) so we may assume that the former can be neglected (it is a generalization of the 

instantaneous approximation in QCD). Let us introduce a vector 


no 


^sinhr;, cosh r/, 0j_ 


(A7) 


where 77 = In 



and make a Sudakov decomposition of the vector < 7 , 


, Po . 

q = A-b rno + q± 

m q 


(AS) 


In a rest frame the projection of the vector q onto po is suppressed as O ( a s ) compared to the other components, for 
this reason in a moving frame we expect that dependence on A should be negligible, i.e. 

V = V(r,q ± ) (A9) 


We define a wave function as 


^(r, <jj_; P) 


d\ x (q, q- P), 















so (IA5I) can be rewritten as 


/ cLtcLq 

V (t - r-,qj_ -p±) ^(t, gj_; P) 


(A10) 


(All) 


and r_ = p_ • uq. Defining A_ = p- ■ po/m q and taking the integral over it in both parts of (IA10I) . we may get 

( B - £ - (r ~ T -’ qi - - p) 
where r = —p ■ p 0 = — sinhppo + cosh 77733 . In coordinate space, the corresponding wave function is 

* (C, rj.) = J ^0j± eKT+ir± ' q± 'lHr, 9 ±; P) 

where the parameter £ is related to ordinary coordinates as 

C, = x o sinh rj — X 3 cosh 77 . 

and has a meaning of ^-coordinate in the J/’F rest frame. The corresponding eigenvalue equation has a form of a 
Schroedinger equation in the rest frame, 


( E + ^U ( C , rj _) = V (C, r±) vh (C, rj_), 

y m q m q J 

and thus guarantees a correct spectrum. For a nonstationary state, we have to replace E —> id ^ 1 where 

(; = x 3 sinh q — Xq cosh 77 , 

so the equation of motion is 


(A12) 


Replacing 


+ — - — ] 'F (C, r±) = V (C, r x ) ’F (C, Pl), 

m q m q I 


= — (cosh r/do + sinh 77 ^ 3 ), 
= — (sinh rjdo + cosh 77 ^ 3 ), 


(A13) 


in the lab-frame, we end up with the second order (w.r.t. time) equation 


-i (cosh,a, + sinh,*) + ^ sin l 1 ’i^ n + cosh ,*) 2 _ £\ „ = y * 

m q m q 


(A14) 


Appendix B: Derivation of (1241251) 


As was discussed in section El in order to find a relation between the rest frame and light front potentials, it is 
necessary to rewrite the gluon momentum in terms of the light-front components, as 


k 2 = ( q-p) =(x-xi) 


m q + P± 




X\ 


■fli 


XX\ 


(x-xi f -q 2 ± + 0 (a 3 s ) . 


(Bl) 


where for the sake of brevity we introduced a notation k = q — p for the gluon momentum. If *F is the light-front wave 
function, and *F (xi,p±) is its Fourier transform over the transverse components, for the interaction term we get 
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U'i = 


J (27T)2 e J 


W 

d 2 p± 

( 2 tt )2 




<B2) 


(*L + /f (x- x,) 2 ) * fe.Px) - 


r d 2 k> sin ( r y k l + ^t^ x ~ x O , 

= 2 J dxi r 2 dr -—- L-\r(T\o ,fcj -' r - L 


y(r)e“ ifc - L ' r - L ^(a;i,r±) = 


r \J k l + ^t( x ~ x O 

= J dxiK (x,xi,r±) 4' {xi,r ±), 
where we introduced a kernel K defined as 


K (x, xi,r±) = 2 J 


k±dk± 

2 n 


dr ■ 


r sin l r ^/kl + ^(x- Xl y 


V(r)J 0 (fcxrj.) 


(B3) 


k l + TrA x ~ *i) 

In the nonrelativistic limit we may approximate * xi ~ 1/4 and simplify the kernel (II j. 'll) to 


K 


n OO pC 

(x,xi,rx)~2 / drrV{r) 

Jo Jo 


k± dk± 

2 tt 


sin ^ r-\Jk\ + 4 m 2 (x — Xi) 2 ^j 
\Jk\+ 4 m 2 (x - xi) 2 


Jo (rx^x) ■ 


(B4) 


The integral over fcxmay be taken introducing a new variable k = yJk\ + 4?n 2 (a; — xi) 2 and using Integral 6.677.1 
from Ryzhik, Gradstein, as 


= 2 

1 

7r 


:(x,xi,rx) « 2 / drrV{r) 

Jo J 2 

r oo 

/ drrV (r) 

Jo 

/*oo 

/ drrV (r) 

J r i 




2m 9 |ai—xi| 

e (r > r±) cos ( 2to 9 I* - 1 \A 2 “’x) 

cos ^2ro g |x — xi| 


sin (rft) Jq ^rx \Jk 2 ~ 4to 2 (x — xi ) 2 ^ = 




(B5) 

(B6) 

(B7) 


Now again change a variable of integration, 2 = 1 /r 2 — r^_, dz = rdr/^r 2 — 7j_, z € (0, 00 ), r = \/r\ + z 2 , so we get 
a result in a very clear and elegant form: 


=> K (x, xi, r±) = — J dzV ^Jr\+ z 2 ^ cos (2m q |x — xi | z), 
(x,rx) = /jxiA'(x,xi,rx)^(xi,xx). 


(B8) 

(B9) 


Notice that since cos(...) is an even function, we may omit an absolute value sign, |x — xi| — > x — x\. Also, using 
symmetry w.r.t. 2 —> — 2 , we can extend the integration region from —00 to + 00 , adding an extra prefactor 1/2, 
and also adding i f_ °° Jzsin {2m |x — Xi| 2), which after integration of a symmetric potential yields zero, so the final 
result for the kernel K may be cast into an equivalent form 


K (x,xi,rx) 



exp (2 im q (x — xi) 2 ). 


(BIO) 
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